Spectral shapes of solid neon 
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We present a Path Integral Monte Carlo calculation of the first three moments of the displacement- 
displacement correlation functions of solid neon at different temperatures for longitudinal and trans- 
verse phonon modes. The Lennard- Jones potential is considered. The relevance of the quantum 
effects on the frequency position of the peak and principally on the line- width of the spectral shape 
is clearly pointed out. The spectrum is reconstructed via a continued fraction expansion; the ap- 
proximations introduced using the effective potential quantum molecular dynamics are discussed. 



Rare gas solids (RGS) are the simplest real systems 
in which we can study lattice vibrations. Argon and 
the heavier RGS can be well approximated as a set of 
harmonic oscillators at lowest temperatures, while the 
classical behaviour is reached before the melting point. 
Quantum corrections on the thermodynamic quantities, 
like kinetic energy and specific heat, can be taken into a< 
count by means of the effective potential approach (EP) 
only for small quantum coupling g, (g < 0.25) defined 
as the ratio between the characteristic frequency and the 
strength of the binding potential. For neon (g = 0.694) 
anharmonic effects are present even for determining the 
ground stateatl. Indeed, precise calculations of kinetic 
energy, to be compared with accurate experiments done 
by Deep Inelastic Neutron Scattering (DINSp, needed a 
rather sophisticated Path Integral Monte Carlo (PIMC) 
computation, being the EP approach inadequate. These 
results show that quantum effects are very important also 
at rather high temperaturesH. 

Information about phonon dynamics is given by spec- 
tral shape, namely the space and time Fourier transform 
of the (symmetrized) displacement-displacement correla- 
tion function: 
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xf is a-th component of the displacement of the i-th atom 
from its equilibrium position. Even though this quantity 
has been investigated since long timea, complete informa- 
tion is not available for various rare gas solids and in par- 
ticular for neon. The perturbative many-body approach 
can give frequency and lifetime of phonons only at low 
temperatures. Classical molecular dynamics (CMD) can 
describe the behaviour of argon and krypton at highest 
temperatures, but it is no more valid for lower tempera- 
tures or stronger coupling. As shown in the following, the 



spectra of neon present significant quantum effects up to 
the melting point and the high quantum coupling pre- 
vents us to use the EP method in the entire temperature 
range. 

We approach the calculation of the spectra of neon as 
given in equation ([[]), by PIMC, evaluating the first three 
even frequency moments, 
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while the odd moments vanish for symmetry reasons. As 
it is well known, this involve the PIMC calculation of 
static correlations obtained by multiple commutators of 
x%(t) = N~ 1/2 Y<iexp(ik- i)xf(t) with the Hamiltonian: 



d n x P , 



dt n 



d n x 



dt n 



(4) 



with the derivatives taken at t = 0. Here we will refer 
only to k = 27r/ao(l, 0, 0) for which one longitudinal and 
two degenerate transverse modes are present, so that we 
omit polarization indexes. 

When the spectra are sufficiently narrow, the normal- 
ized second moment, Sik and the irreducible part of 
fourth moment, 62k- 
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can be directlK related with the peak position and width 
of the spectraQ. 

A reconstruction of the spectra can be done by the 
continued fraction expansion of the Laplace transform of 
the normalized correlation function, 
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with 



S s (k, u) = 5RS (k, z = iuj) C(k, t = 0). 
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The continued fraction expansion can be stopped at the 
third stage with a suitable termination S2(k, z) 
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In this paper, we shall present some of these spectra, 
showing the validity of the approach. Neutron scattering 
data are not available up to now and we suggest and dis- 
cuss here the possibility to perform such an experiment. 

We have considered samples of solid neon, with 256 
atoms interacting through a (12-6) Lennard- Jones pair- 
wise potential, 



V(r) = 4e 



(9) 



and periodic boundary conditions. The dynamic interac- 
tion is limited to the 12 nearest neighbours, while for the 
outer shells the static approximation is used. For every 
temperature we resorted to three different Trotter num- 
bers: P = 8, 16, 24. For each of them 16 simulation runs 
of 100 000 steps per particle were performed, plus 20 000 
steps per particle for initial thermalization. The den- 
sity was adjusted in order to get a practically vanishing 
pressure (the pressure is always less than 15 atm). The 
parameters of the Lennard- JUj>nes potential are taken as 
e = 36.68 K and a = 2.787 All The melting temperature 
of neon at zero pressure is 24.5 K. 

Detailed explanations of the total .-procedure and re- 
sults for the first three even moments^, at different tem- 
peratures for longitudinal and transverse modes, will 
be presented in an extended paper. Here we want to 
point out that when the order of the frequency moment 
increases, more complicated static correlations are in- 
volved. Moreover accurate Trotter extrapolations are in 
order and finite-size effects are more and more impc 
tant. We have used the procedure introduced by us 
by which we correct the raw PIMC data, subtracting the 
exact contributions of the harmonic part for finite P and 
N, and adding both the exact harmonic results for infi- 
nite P and N, In this way, an accuracy of 0.2% is reached 
for the zeroth moment, which rises to 1% for the fourth 
moment. This corresponds to a maximum uncertainty 
of 4% for i$2k- Classical simulations were also done for 
comparison, using both the classical and the effective po- 
tential. 

The necessity to account for quantum effects by an 
"exact" method like PIMC has been ascertained, espe- 
cially for higher-order moments. The fourth moment ex- 
hibits strong quantum effects at any realistic temperature 
which cannot be approached by the EP method. In tables 
I and H we report the first three even moments together 
with the analogous ones obtained by classical and EP 
Monte Carlo calculations. The fourth moment is related 
to the width of the spectra (phonon lifetime) through 
the quantity 62k- For narrow spectra ((52k <C #ik), the 
phonon frequencies are cjk ~ V^ik while the phonon life- 
times are 77k ~ \/^2k/2- These quantities, shown in table 
III, can be probed by inelastic neutron scattering. 




FIG. 1. Trasverse projection of the (normalized) space 
and time Fourier transform of the symmetrized displace- 
ment-displacement correlation function, 5s (k, ui)/C(k, t — 0) . 
at T = 10 K p = 1.494 g cm~ 3 for k = 27r/a (l, 0, 0) and 
zero pressure. The dashed line refers to CMD results, the 
long-dashed one to EPMD results, and the solid line to re- 
sults obtained via continued fraction expansion and PIMC 
evaluation of frequency moments. 




(0 (meV) 

FIG. 2. The same as in figure |l| for the longitudinal projec- 
tion of the (normalized) space and time Fourier transform of 
the symmetrized displacement-displacement correlation func- 
tion, S s (k,w)/C(k,t = 0), at T = 10 K, p = 1.494 g cm~ 3 
for k = 27r/a (l,0,0) 

Finally, we have calculated the spectra with the avail- 
able <5ik, <^2k by means of the continued fraction (||) and 
with a suitable "Gaussian termination"ll3: 32(k, t) — 
exp (— rt 2 ) . The parameter T, is determined by the in- 
sight on the corresponding spectra obtained by CMLliil. 

As an example, some spectra are shown in figures |l| 
and ^|. In these figures we report also the similar spectra, 
obtained by EP molecular dynamics (EPMD)E3. 

We therefore can conclude that: 

i) Quantum effects in solid neon are relevant at all tem- 



2 



peratures so that neon cannot be approached by 
classical models. 

ii) The evaluation of the spectral width requires a par- 
ticular care and a fully quantum treatment of the 
fourth moment. EPMD, as observed int2l, can cor- 
rectly give the peak position only because it re- 
produces just the short-time behaviour and conse- 
quently the second moment; using this method the 
phonon damping is therefore calculated by consid- 
ering classical processes only. This is not sufficient 
for giving a good description of the spectra as on 
the other hand our approach is expected to do. 

We conclude suggesting new accurate neutron scattering 
experiments in order to investigate quantum effects in 
solid neon. In particular our results on zero moment 
can be tested measuring the integrated intensity, while 
the afore-mentioned features of the line shapes can be 
directly compared with experimental spectra. 

We would like to thank Prof. S. W. Lovesey for useful 
discussions. One of us (V. T.) wants also to thank him 
for hospitality in the very stimulating atmosphere of the 
Rutherford Appleton Laboratory. 
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TABLE I. Moments of S s (k,u) for T = 10 K at vanishing 
pressure, for the wave vector k = 2n/ao(l, 0, 0). They are 
expressed in reduced units, b — 2.787 A and ujq = 0.289747 
meV. Classical second moments are exact since they are equal 
to the reduced temperature. 
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98.13 ±0.05 


168.5 ±0.1 


314. ±3. 
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M /6 2 (xhn 4 ) 


3.955 ±0.002 


20.732 ± 0.008 


19.02 ±0.05 


M 2 /(b 2 Lul) 


0.2726281 


0.9902 ± 0.0003 


1.141 ±0.004 


M 4 /(b 2 ul) 


205.0 ± 0.2 


467.3 ± 0.3 


849. ± 6. 



TABLE II. Moments of 5 s (k, uj) for T = 20 K at vanishing 
pressure, for the wave vector k = 27r/ao(l, 0, 0). They are 
expressed in reduced units, b — 2.787 A and ujq = 0.289747 
meV. Classical second moments are exact since they are equal 
to the reduced temperature. 







Transverse 






classical 


EP 


PIMC 


M /6 2 (xhn 4 ) 


20.92 ±0.02 


36.50 ± 0.05 


34.7 ±0.3 


M 2 /(6 2 ^ 2 ) 


0.545256 


0.8387 ± 0.0006 


0.870 ± 0.002 


M 4 /(6 2 ^) 


182.3 ±0.3 


218.6 ±0.4 


347. ± 3. 
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9.67 ±0.01 


22.47 ±0.02 


21.3 ±0.1 


M 2 /(6 2 ^ 2 ) 


0.545256 


1.1022 ±0.0006 


1.157 ±0.006 


M 4 /(6 2 w 4 ) 


380.3 ±0.5 


560.1 ±0.3 


872. ± 8. 



TABLE III. v 7 ^! and Vw/2 roughly represent the peak 
position and the phonon life time, respectively. They are 
evaluated via PIMC simulations and they are expressed in 
meV. 
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10 


4.83 ±0.02 1.58 ±0.1 


7.01 ±0.02 1.78 ±0.06 


20 


4.57 ±0.02 1.74 ±0.1 


6.76 ±0.02 2.07 ±0.06 
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